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Finding the factors of an integer can be achieved by various experimental techniques, based on an 
algorithm developed by Schleich et al., which uses specific properties of Gaufi sums. Experimental 
limitations usually require truncation of these series, but if the truncation parameter is too small, it 
is no longer possible to distinguish between factors and so-called "ghost" factors. Here, we discuss 
two techniques for distinguishing between true factors and ghost factors while keeping the number 
of terms in the sum constant or only slowly increasing. We experimentally test these modified 
algorithms in a nuclear spin system, using NMR. 
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Introduction — Factorization of large numbers is a 
computationally hard problem: the computational re- 
sources required to accomplish this task increase expo- 
nentially with the number of digits [1] for all algorithms 
discovered until 1994. Then Peter Shor developed an al- 
gorithm that can solve the task in polynomial time. This 
algorithm requires a computational device that operates 
according to the laws of quantum mechanics, storing in- 
formation in quantum states and performing logical oper- 
ations as unitary evolutions under suitable Hamiltonians 
2 . Experimental implementations of Shor's factoriza- 
tion algorithm were demonstrated first with nuclear spins 
as qubits [5] , and recently with photonic qubits [U [5] . 

More recently, another factorization algorithm was 
proposed by Schleich and co-workers [5J [HI IH] which 
uses properties of GauB sums. A complete normalized 
quadratic GauB sum is defined by 
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where N is the integer to be factorized and I is the trial 
factor. If Ms a factor of N, i.e., N/l is an integer, 
the resulting sum is 

I^^OI = 1- In all othe r cases, 

\A l N \l)\ < 1. 

The number of terms that has., to be evaluated for 
the complete GauB sum of Eq. (HI grows as ' = 

\^/N{-\/N — 1) oc N. A factorization algorithm on the 
basis of Eq. (|T|) is thus computationally very expen- 
sive. However, in most CELS6S, cl complete evaluation is 
not necessary. Recent experimental implementations us- 
ing NMR [Till E] , cold atoms [T^] and ultra-short laser 
pulses [13 have successfully demonstrated that it is usu- 
ally possible to truncate the sums after a relatively small 
number of terms. We write the corresponding truncated 
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with a constant truncation parameter M for each argu- 
ment /, instead of the upper limit I — 1 in the complete 
Gaufi sum of Eq. Accordingly, only M\fN terms 

have to be added, greatly improves the efficiency and 
precision of the experiments. However, the truncation of 
the Gaufi sum weakens the discrimination of the factors 
from nonfactors, resulting in the appearance of "ghost" 
factors, whose Gaufi sums are close to unity. The re- 
quirement of suppressing these "ghost" factors thus sets 
a lower limit on the choice of the truncation parameter 
M. 

The choice of an optimal truncation parameter M was 
discussed in a recent paper by Stefanak et al. |14j . 
They found an upper bound on the truncation parame- 
ter M ~ \/N, which represents a sufficient and necessary 
condition for the success of the Gaufi sum factorization 
scheme. Here, we experimentally explore this issue, using 
liquid-state NMR for the evaluation of the Gaufi sums. 
Furthermore, we demonstrate additional possible ways 
of reducing the truncation parameter M, while keeping 
excellent contrast between factors and non-factors. 

Effect of truncation — The choice of the truncation pa- 
rameter M plays a crucial role in the success of the Gaufi 
sum factorization scheme. In the previous experiments, 
the visibility of the resulting factorization interference 
pattern was high enough for successful factorization of, 
e.g., the eight-digit number N — 52882363 by logarith- 
mically choosing the truncation parameter M = InN 
PHI [IT]. However, in some cases, we also observed 
some "ghost" factors. Like in Ref. [14j. we define those 
trial factors for which the absolute value of the trun- 
cated Gaufi sum is larger than the threshold value 1 / \/2 
as "ghost" factors. 

As Stefanak et al. showed [13] . this sum behavior for 
different trial factors is best analyzed by considering the 
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fractional part of 2N/1, 
e(N,l) 
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with |e| < 1. Here 2k is the closest even integer to 2N/1. 
Since exp(i2nm 2 k) = 1, the GauB sum ^ can be rewrit- 
ten as 



A%{1) = s M (e) 
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where sm(c) is the normalized curlicue function, which 
has the property: 

!1, e = for all factors, 
-j=, e = 0.5 for threshold non-factors{5) 
0, e — > 1 for typical non-factors. 

Here three different classes for the trial factors are de- 
fined [3]. For the class of the ghost factors, the curlicue 
function depends on the truncation parameter M: 



1, for a small M, 

0, for a very large M. 



(6) 



Ghost factors occur when e is very close to zero and the 
GauB sum is truncated after too few terms. Figure [l] 
illustrates this behavior: for a given e, the number of 
terms M needed to suppress the value of %(e) below 
the threshold l/y/2 is « 1/ 
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FIG. 1: (color online). Absolute values |sm(c)| of the nor- 
malized curlicue function vs. the truncation parameter M for 
different values of e. 

Fig. [2] illustrates how ghost factors occur for small 
e: the three parts of the figure show the distribution 
of the different terms of a GauB sum in the complex 
plane for (a) M = 20, (b) 200, and (c) 1000 terms for 
e = 4x 10~ 5 : only in the third case are the phases of the 
terms sufficiently well distributed that the sum (shown 
as the blue line) approaches zero. 

Randomized procedure — Suppressing all ghost factors 
of truncated GauB sums below the threshold \j\/2, re- 
quires M ~ [15]. This may still be too large for ex- 
perimental factorization of large integers, e.g., M ~ 1000 




FIG. 2: (color online). Distribution of the phases <f> m — Tvm 2 £ 
for a fractional part e(N, 1) = 4 X 10~ 5 for different truncation 
parameters M: M = 20, M = 200 and M = 1000. The 
individual terms of the GauB sums are represented by red 
dots. The arrows from the origin (0, 0) represent the resulting 
truncated Gaufi sums of Eq. pj, whose absolute values are, 
respectively, |s A /=2o| ~ l,|sjwr=20o| ~ 0.3155 and |sm=iooo| ~ 
0.0770. In (c), the black points on the circle are the random 
phases created by a randomized procedure: 10 values of m 
were randomly chosen from the range [0, M m ax]- Their sum 
is ~ 0.0023, very close to the origin. 



for the factorization of a 12-digit integer requires ex- 
treme precision in the experimental implementation and 
may lead to excessive decoherence. However, the num- 
ber of terms in the sum can be reduced significantly be- 
low the v^/V threshold by a suitable choice of the terms 
that are evaluated: instead of evaluating all terms with 
< m < M, we randomly choose a relatively small frac- 
tion of the terms with m 6 [0, M ~ \/N]. 

As a test of this procedure, we select a "difficult" case, 
where the conventional procedure requires a large trun- 
cation parameter. Such cases occur, e.g., for numbers 
that are products of neighboring primes. We chose as an 
experimental example the product of the 100000 th and 
prime, 



100001 th 



N = 1689259081189 = 1299709 x 1299721. 

In the experiments, we implement the procedure by 
adding nuclear spin coherence, using the ^ (7 = |) nu- 
clear spins of water, diluted in D 2 0. The nuclear system 
was contolled by suitable radio-frequency (rf) magnetic 
fields. In the rotating frame, an rf pulse with duration r, 
amplitude u and phase cj> m generates the unitary opera- 
tor 



U m = exp(-i9(I x cos <j> m + I y sin <j> m )). 



(7) 



The lower index m indicates that the Hamiltonian is 
specific for each term in the GauB sum. The flip an- 
gle 9 — lot represents the absolute value and </> m the 
phase of a complex number in the series, corresponding 
to <j> m (l) = 2i:m 2 ^Y in the GauB sum. The sum was real- 
ized in the experiment by applying a sequence of M+l rf 
pulses with small flip-angle to the spins, with the phase 
of each pulse equal to that of the corresponding term of 
the GauB sum. A short delay (5/js) was inserted between 
the pulses. The combined effect of these pulses can be 
described by the propagator U(l) = Um • • • U$. 

In the limit of small flip angles, M9 <C 1, the oper- 
ators in the exponent approximately commute and the 
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propagator can be approximated by 
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If i is a factor of N, then 4> m (l) = 2kir with /c integer and 
all M+l pulses have the same phase (0 = 0). In this case, 
the combined effect of the pulses is U{1) ss e -^(M+i)i x _ 
If it is applied to the thermal equilibrium state, it cre- 
ates transverse magnetization I y , with an amplitude 
oc 9(M + 1). If Z is not a factor of N, the individual 
signals interfere destructively and the resulting trans- 
verse magnetization is close to zero. For each experiment, 
the transverse magnetization generated was recorded as 
a free induction decay (FID). 

The experiments were carried out on a 500 MHz Bruker 
Avance 11+ NMR spectrometer. Using the standard 
truncated GauB sum Aj$(l) of Eq. ^ with M = 19, 
we obtained experimental results that were indistinguish- 
able from the maximal value of 1 for all trial factors close 
to the real factors [see upper trace in Fig. [3]. However, 
if we use the randomly selected m-values, as few as 10 
terms are sufficient to suppress all the non-factors well 
below the threshold of 1/ \/2, as shown in Fig. [3] (blue 
spectra and red dots), while the real factors 1299709 and 
1299721 always yielded values close to 1. 

As the second example, we chose to factorize a 17-digit 
integer 

N = 32193216510801043 = 179424673 x 179424691. 

Randomly choosing 10 values of m from the range 
[0,M max = 5000], we experimentally evaluated the par- 
tial GauB sums for the trial factors I between 179424663 
and 179424701. The results, shown in Fig. [4] clearly 
show that the factors 179424673 and 129424691 are found 
and no ghost factors appear. 

Exponential sums — Exponential sums can be defined 

as 
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In terms of the fractional part e(N, L) of 2N/1, they are 
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The standard case is recovered for n — 2. These higher- 
order exponential sums can be used for factorization ex- 
actly as the second-order function: again, the factors 
generate constant phases for all terms and thus the maxi- 
mum of the sum, while non-factors ideally generate sums 
much smaller than unity. These higher-order exponential 
sums can provide higher contrast between factors and 
non-factors, even for small truncation parameters M. As 
shown in Fig. [5j the higher the order n, the smaller the 
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FIG. 3: (color online). Factorization of N = 1689259081189 
with tmin ~ 1.693xl0~ 5 . The upper trace shows the standard 
truncated Gaufi sum with the truncation parameter M = 19, 
the middle trace shows the result of the Monte Carlo proce- 
dure with 10 randomly chosen values of m from the range 
[0,M ma x = 1000], and the lower trace shows the 5 th -order 
truncated exponential sum A5f} with the truncation param- 
eter M = 10. While all trial factors masked as true factors 
in the standard truncated procedure (upper trace), the true 
factors are easily found in the Monte Carlo and exponential 
sum procedure. 



upper bound of M required to suppress the value of | 
below the threshold. 

Numerical analysis show that the truncation parame- 
ter M required to push all non-factors below the thresh- 
old l/\/2 scales with the order n of the exponential 
sum and the size of the number N to be factorized as 
M ~ 2 \/N [15]. Therefore, to factorize a 12-digit inte- 
ger, the required value of M decreases from 10 3 to 10 
if we use the 6 th order function instead of the quadratic 
truncated GauB sum. However the advantage takes the 
price of the smaller gap between factors and threshold 
non- factors [TH]. The authors also proposed an NMR re- 
alization of exponential sums |16j . 

We experimentally tested the performance of the 
higher-order truncated exponential sums A'nM , using the 
same procedure as for the n = 2 case. The lower trace 
of Fig. [3] clearly shows that this procedure provides an 
excellent contrast between factors and non-factors, even 
for a relatively small number of terms M — 10. 

Conclusion — GauB sums [T71 [TH1 [25] are ubiquitous 
in number theory and found many applications, such as 
Plancherel's theorem on finite groups [20], the Talbot ef- 
fect of classical optics [2T], fractional revivals [251 l2"3"] . 
quantum carpets [M] and Josephson junctions j^S]. Re- 
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FIG. 4: (color online). Factorization of N = 
32193216510801043 by the randomized procedure with 10 ran- 
dom phases in [0, M max = 5000]. 
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cently, Gaufi sums were also used for factorization, which 
is related to the proposal of Clauser and Dowling [26] to 
factor an integer using a familiar Youngs N-slit classical 
interferometer . 



In this paper, we presented an experimental investiga- 
tion on the GauB sum factorization scheme for large num- 
bers, where ghost factors often appear when the trunca- 
tion parameter M is relatively small (e.g., M ~ 15 — 20). 
In these cases, the truncation parameter M must be in- 
creased to relatively large numbers, which is undesirable 
for experimental implementations. To circumvent this 
increase in the required number of terms, we have in- 
troduced a Monte Carlo procedure, where the required 
number of terms remains roughly constant, and have 
used higher-order truncated Gaufi sums, whose scaling 
behavior is much more benign than for the second order 
function. While we have used a nuclear spin system for 
the experimental implementation, it should be straight- 
forward to apply this scheme to other (quantum or clas- 
sical) systems. 
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FIG. 5: (color online). Comparison the different orders of 
Gaufi sums as a function of the truncation parameter M for 
a fractional part e = 1.0 x 10 -6 . 



Note added: During the preparation of this paper, we 
became aware of closely related work J27J [28] . 



[1] D. E. Knuth, The Art of Computer Programming Vol. 

2, Seminumerical Algorithms (Addison- Wesley, Reading, 

Massachusetts, 1998). 
[2] P. Shor, in Proceedings of the 35th Annual Symposium 

on Foundations of Computer Science (IEEE Computer 

Society Press, New York, 1994), p. 124. 
[3] L. M. K. Vandersypen, M. Steffen, G. Breyta, C. S. Yan- 

noni, M. H. Sherwood, and I. L. Chuang, Nature 414, 

883 (2001). 

[4] C.-Y. Lu, D. E. Browne, T. Yang, and J.-W. Pan, Phys. 
Rev. Lett. 99, 250504 (2007). 

[5] B. P. Lanyon, T. J. Weinhold, N. K. Langford, M. Bar- 
bieri, D. F. James, A. Gilchrist, and A. G. White, Phys. 
Rev. Lett. 99, 250505 (2007). 

[6] W.Merkel, O. Crasser, F. Haug, E. Lutz, H. Mack, 
M. Freyberger, W. P. Schleich, I. Averbukh, M. Bienert, 
B. Girard, et al., Int. J. Mod. Phys. B 20, 1893 (2006). 

[7] W.Merkel, I. S. Averbukh, B. Girard, G. G. Paulus, and 
W. P. Schleich, Fortsch. Phys. 54, 856 (2006). 

[8] W. Merkel, W. P. Schleich, I. S. Averbukh, and B. Gi- 
rard, unpublished. 



[9] W. Merkel, W. P. Schleich, G. G. Paulus, and B. Girard, 
unpublised. 

[10] T. S. Mahesh, N. Rajendran, X. Peng, and D. Suter, 

Phys. Rev. A 75, 062303 (2007). 
[11] M. Mehring, K. Miiller, I. S. Averbukh, W. Merkel, and 

W. P. Schleich, Phys. Rev. Lett. 98, 120502 (2007). 
[12] M. Gilowski, T. Wendrich, T. Miiller, C. Jentsch, 

E. M. Rasel, W. Ertmer, and W. P. Schleich, 

arXiv:0709.1424vl, accepted in Phys. Rev. Lett. (2007). 
[13] D. Bigourd, B. Chatel, B. Girard, and W. P. Schleich, 

arXiv:0709.1906vl, accepted in Phys. Rev. Lett. (2007). 
[14] M. Stefanak, W. Merkel, W. P. Schleich, D. Haase, and 

H. Maier, New J. Phys. 9, 370 (2007). 
[15] M. Stefanak, D. Haase, W. Merkel, M. S. Zubairy, and 

W. P. Schleich (2008), submitted to J. Phys. A. 
[16] M. Stefanak, W. Merkel, M. Mehring, and W. P. Schleich 

(2008), submitted to J. Mod. Phys. 
[17] S. Lang, Algebraic Number Theory (New York: Addison 

Wesley, 1970). 

[18] H. Maier and W. P. Schleich, Prime Numbers 101: 
A Primer on Number Theory (New York:Wiley-VCH, 



5 

2007). 

[19] H. Davenport, Multiplicative Number Theory (New York: 

Springer, 1980). 
[20] K. Yosida, Functional Analysis (Springer Verlag, 1968). 
[21] H. F. Talbot, Phil. Mag. 9, 401 (1836). 
[22] C. Leichtle, I. S. Avcrbukh, and W. P. Schleich, Phys. 

Rev. A 54, 5299 (1996). 
[23] C. Leichtle, I. S. Averbukh, and W. P. Schleich, Phys. 

Rev. Lett. 77, 3999 (1996). 
[24] M. V. Berry, I. Marzoli, and W. P. Schleich, Phys. World 

14, 39 (2001). 

[25] J. Oppenlander, C. Haussler, and N. Schopohl, Phys. 

Rev. B 63, 024511 (2000). 
[26] J. F. Clauser and J. P. Dowling, Phys. Rev. A 53, 4587 

(1996). 

[27] E. Rasel, 38th Winter Colloquium on the Physics of 

Quantum Electronics, Snowbird, US (Jan 2008). 
[28] B. Chatel, 38th Winter Colloquium on the Physics of 

Quantum Electronics, Snowbird, US (Jan 2008). 



